Simplest nonequilibrium phase transition into an absorbing state 
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We study in further detail particle models displaying a boundary-induced absorbing state phase 
transition [Phys. Rev. E. 65, 046104 (2002) and Phys. Rev. Lett. 100, 165701 (2008)] . These are 
one-dimensional systems consisting of a single site (the boundary) where creation and annihilation 
of particles occur, and a bulk where particles move diffusively. We study different versions of these 
models, and confirm that, except for one exactly solvable bosonic variant exhibiting a discontinuous 
transition and trivial exponents, all the others display non-trivial behavior, with critical exponents 
differing from their mean-field values, representing a universality class. Finally, the relation of these 
systems with a (0 + 1) -dimensional non-Markovian process is discussed. 

PACS numbers: 64.60.Ht, 68.35.Rh, 64.70.-p 



I. INTRODUCTION 

Phase transitions occurring in the bulk, but driven by 
specific conditions at its boundaries, are called boundary- 
induced phase transitions [l| . Examples include diffusive 
transport @, [|| and traffic flow Q models. A simple ex- 
ample for this is provided by the one-dimensional totally 
asymmetric simple exclusion process Q, where particles 
enter the system at the left boundary, jump to the right 
in the bulk, and exit at the right boundary. Depending 
on the entering and exiting rate values, the system ex- 
hibits qualitatively different phenomenologies (maximal 
current, large current and low density, or small current 
and high density), with straightforward applications to 
traffic flow problems. 

In the present work, we are interested in boundary- 
induced phase transitions in systems with absorbing 
states. An absorbing state is a dynamical trap which 
can be accessed but cannot be left d, 0, [![. Systems 
with absorbing phase transitions are controlled by a pa- 
rameter, depending on which the system either enters the 
absorbing state with certainty or survives in a stationary 
fluctuating/active state. The most prominent family of 
phase transitions into an absorbing state is the very ro- 
bust direct percolation (DP) universality class. A recent 
breakthrough has been the experimental observation of 
DP critical behavior for the first time @. 

A paradigmatic model in the DP universality class is 
the contact process (CP) It can be viewed as a sim- 
ple model for the propagation of a disease where sick in- 
dividuals can infect healthy neighbors or become healthy 
spontaneously. More precisely, in the CP in d spatial 
dimensions, a particle (infected individual) can be cre- 
ated at an "empty" site with a rate Xn/2d, where n is 
the number of nearest neighbors occupied by a particle, 
and an occupied site can become empty at rate 1. The 
empty configuration is an absorbing state. For A larger 



than a certain critical threshold, A c , the process is able 
to sustain (in an infinite lattice) a non-vanishing den- 
sity of particles, while for A < A c the dynamics ends up, 
ineluctably, in the absorbing state. 

As continuous phase transitions involve long-range cor- 
relations, boundary effects may play an important role. 
In the context of absorbing phase transitions, previous 
studies focused primaril y o n DP confined to parabo- 
las [111, IT3 | , ac tive walls |13j. as well as absorbing walls 
and edges [3, 0j| . Although such boundaries influence 
the dynamics deep into the bulk, the universality class 
of the bulk transition is not inherently changed, rather it 
is extended by an additional independent exponent de- 
scribing the order parameter near the boundary. There- 
fore, the question arises whether it is possible to find 
boundary-induced absorbing phase transitions, absent in 
the corresponding systems without boundaries, consti- 
tuting independent universality classes. 

In this paper we present a detailed discussion of two 
slightly different models introduced in Ref. [l6[ and [ItJ , 
respectively. Both of them exhibit a boundary-induced 
nonequilibrium phase transition into an absorbing state. 
These are one-dimensional particle systems consisting of 
a single site (and, at most, its nearest neighbor), where 
creation and annihilation of particles occur, and a bulk, 
where particles move diffusively. While in the first ref- 
erence [l6|], the dynamics at the boundary is a contact 
process, in the second one [lj} particles at the origin an- 
nihilate only pairwise (A + A — > 0). We study different 
versions of these models to compare them and scruti- 
nize the relevance of relaxing the fermionic constraint 
(i.e. occupation number not restricted to or 1) both at 
the bulk and at the boundary. Our study includes mean 
field approximations, numerical analysis, some field the- 
oretical arguments, as well as the relation with a (0+1)- 
dimensional non-Markovian model [i~8j . 

The paper is organized as follows: In the next section 
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FIG. 1: (Color online) Model on a semi-infinite lattice with 
symmetric diffusion in the bulk and special dynamical rules 
at the left boundary (see main text). 

we define the first model and present numerical results. 
In Sec. IIIII we discuss various types of mean field approx- 
imations and show that this model has indeed a non- 
trivial behavior. Sec. IIVI is concerned with two bosonic 
versions of the first model and the study of models with 
pair annihilation at the boundary. One of the bosonic 
versions is solved exactly and it is shown to have triv- 
ial critical behavior. Instead, the other bosonic version 
and models with pair annihilation are shown to share 
the same critical behavior as the first model. In Sec. [V] 
we discuss the relation with a (0+l)-dimensional non- 
Mar kovian model [l8| and, finally, we present our main 
conclusions. 
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FIG. 2: (Color online) Density of particles at the leftmost site 
po (left) and the average total number of particles TV (right) 
as functions of time for different values of A, below, above, 
and at criticality. At the critical point, po(t) decays as t -1 ^ 2 
while and (N(t)) is essentially constant. 

at p c — 0.74435(15) has been reported to separate the 
absorbing phase, in which the total number of particles 
vanishes, from another with indefinitely sustained activ- 
ity [l6[. 



II. BASIC MODEL DEFINITION AND 
SIMULATIONS 

A. Definition of the model 

The model presented in [l6| is defined on a one- 
dimensional semi-infinite discrete lattice where each site 
is either occupied by a particle (si — 1) or empty (s, = 0). 
All lattice sites have two neighbors, except for the bound- 
ary [i = 0) with a single one. The dynamics is a com- 
bination of an unbiased random walk in the bulk and a 
contact process-like dynamics at the left boundary. It is 
implemented as follows: 

(a) A particle is randomly selected. 

(b) If it is located at the leftmost site, it generates an- 
other particle at site 1 with probability p, provided 
that it is empty (si = 0), or it dies (so = 0) with 
probability 1 — p. 

(c) Particles in the bulk perform a symmetric exclusion 
process, moving to any of their two neighbors with 
equal probability, provided that the destination site 
is empty (otherwise nothing happens). 

Starting with a single particle at the leftmost site in an 
otherwise absorbing (i.e. empty) configuration, the pro- 
cess evolves as follows: the initial particle at site ei- 
ther dies or generates another particle at the neighbor- 
ing site 1. This last performs a random walk in the bulk 
until, eventually, it returns to the origin to create an- 
other offspring or disappear. A critical point, located 



B. Order parameters 

A possible order parameter for this model is the aver- 
age density of particles at the leftmost site: 

Po = («o) , (1) 

where () stands for ensemble averages. This quantity is 
plotted in Fig. [5] as a function of A := p — p c . At the 
critical point, po(t) decays algebraically in time, as: 

Po (t) ~ r a (2) 

with an exponent a = 0.50(1), compatible with a rational 
value a — 1/2. 

Another possibility is to choose as an order parameter 
the average total number of particles, (N(t)), which, as 
shown in Fig. [21 goes to zero for p < p c and increases 
steadily for p > p c (actually, it is limited only by the 
system size). At criticality, (N(t)) is found to be constant 
in the large time limit. 

In the usual scaling picture of absorbing phase transi- 
tions, the critical exponent (3 is related to the probability 
that a given site belongs to an infinite cluster generated 
from a fully occupied lattice at t = — oo. This quantity 
tends to zero as the control parameter approaches the 
critical value from above. Similarly, the exponent /?' is 
related to the probability that a localized seed generates 
an infinite cluster extending to t = +oo. Therefore, in 
the supercritical phase (A > 0), the averaged activity of 
the site at the origin for t — > oo measured in seed simula- 
tions averaging over all runs, scales as p s ~ A^ + ^ , where 
the superscript 's' stands for 'stationary'. At criticality, 
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FIG. 3: (Color online) Left panel: Density of particles at the 
leftmost site, at the time when it reaches its minimum value, 
as a function of the distance from criticality A. This gives 
the exponent (3 = 0.68(5). Right panel: The same quantity, 
at criticality and in the stationary state, as a function of the 
external field h, giving S^ 1 = 0.29(5), compatible with the 
conjectured value 1/3. 

this function is expected to decay as p(t) ~ t~(P+P )f v \\, 
where vn is the correlation time exponent. Moreover, in 
the DP class a special time reversal symmetry implies 
that (3 = 13' @. 

As shown in [l7j , time reversal symmetry also holds in 
the present type of models. This implies that, in super- 
critical seed simulations, the density of active sites at the 
boundary is expected to saturate as: 

Po~A 2 ^, (3) 

while, at criticality: 

po{t)~t-W"\ (4) 
implying that a in Eq. (J2]) is 

a = 2p/v\\ . (5) 
Assuming that a — 1/2, then (3/v\\ = 1/4. 



Note that this can be true only if the density po(t) in 
seed simulations obeys the scaling relation: 

p (t) = A 2 ? R(tA^) (7) 

i.e. if it is possible to collapse the data by plotting 
PqA~ 213 versus tA 4 @. Indeed, this will be shown to be 
the case in Sec. [V] for a 0-dimensional non-Markovian 
process argued to be in the same universality class. 

Relying on this observation, one can determine the 
value of the exponent (3 by measuring the density 
at the minimum, which is reached much earlier than the 
stationary state. In Fig. [3] we plot p m as a function of A, 
inferring (3 = 0.68(5). 



D. External field 

In ordinary directed percolation, an external field, con- 
jugate to the order parameter, can be implemented by 
creating active sites at some constant rate h, thereby de- 
stroying the absorbing nature of the empty configuration. 
At criticality, the external field is known to drive a tri- 
dimensional DP process towards a stationary state with 
p s ru h 1/Sh where S^ 1 = /3/(v\\ + du± - f3'), and u± is the 
correlation length critical exponent. 

In the present model, the external field, conjugate to 
the order parameter po, corresponds to spontaneous cre- 
ation of activity at the leftmost site at rate h. The above 
hyperscaling relation for Sh is thus expected to be fulfilled 
by taking d = 0: 

Po-h 1 ^. (8) 

with 

s- 1 = /?/h - n (9) 

From this expression, exploiting the fact that f3 = (3' and 
using Eq.([5]) as well as the conjectured rational value a — 
1/2, a prediction S^ 1 = 1/3 is obtained. Our numerical 
estimate, S^ 1 = 0.29(5) (see Fig. [3]) is compatible with 
this result. 



C. Stationary properties 

In numerical simulations in the active phase, it takes a 
very long time, specially for small values of A, to reach 
the steady state. Moreover, we observed the unusual 
fact that, for A > 0, the density po goes through a mini- 
mum before reaching the stationary state (see Fig. [2] and 
also [l9j . where similar non-monotonous curves were re- 
ported). However, it turns out that the value p™ at the 
minimum and the saturation value p s Q differ by a con- 
stant factor, entailing that both quantities scale in the 
same way, i.e.: 

p"o l - A 2/3 . (6) 



E. Survival probability 

The survival probability P s (t) is defined as the fraction 
of runs that, starting with a single seed at the boundary, 
survive at least until time t. At criticality, this quantity 
is expected to decay algebraically: 

Ps{t)~t- 5 , (10) 

with the so-called survival exponent 5, while in the super- 
critical regime it saturates in the long time limit. Since 
P s (oo) coincides with the probability for a seed to gener- 
ate an infinite cluster, the saturation value of the survival 
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FIG. 4: (Color online) Survival probability P s (t) as function 
of time below, above, and at criticality. At the critical point, 
it decays with the exponent S = 0.15(2), different from 
and in agreement with the conjectured value 1/6. 



probability as a function of the distance from criticality 
gives the exponent (3'. As in DP, one expects P s (t) to de- 
cay in time with an exponent 6 — (3' /v\\ = 1/4. However, 
as shown in Fig. [4j one finds a much smaller exponent 
8 = 0.15(2). Therefore, the usual relation 8 = does 
not hold. We also observed that it is not possible to col- 
lapse different curves of P s (t) for different values of A, 
i.e. the survival probability seems to exhibit an anoma- 
lous type of scaling behavior. We expect that off-critical 
simulations of the survival probability give the exponent 
(3' but the simulation times needed to reach steady state 
are prohibitively long. 

An explanation for the value 8 = 0.15(2), differing from 
/3/f\\, is given in the following subsection. 



F. Time reversal symmetry 

In ordinary bond DP, the statistical weight of a con- 
figuration of percolating paths does not depend on the 
direction of time. More specifically, the probability to 
find an open path from at least one site at time t = to 
a particular site at time t coincides with the probability 
to find an open path from a particular site at time t = 
to at least one site at time t. This implies that, in bond 
DP, i) the density p(t) in simulations with fully occupied 
initial state and ii) the survival probability P s (t) in seed 
simulations coincide; hence (3 = (3' . In other realizations 
of DP (e.g. site DP), this time reversal symmetry is not 
exact but only asymptotically realized. 

Applying the same arguments to the present model, 
the survival probability P s (t) in seed simulations should 
scale in the same way as the density of active sites at the 
boundary po [t) in a process starting with a fully occupied 
lattice in the bulk. A numerical test, which approximates 




FIG. 5: (Color online) Left: Data collapse of the rescaled 
profiles of the particle density at criticality for to = 
64, 128, . . . , 8192 (blue) compared to a Gaussian distribution 
(red). Inset: The same data collapse in a double-logarithmic 
representation. Right: Density of particles (blue) and pairs 
(green) at to = 10 6 , showing the presence of correlations 
which decay in space as x~ 1 ' 2 , indicating that /3/u± = 1/2. 



such a situation, confirms this conjecture, i.e. one has 
Po(t) ~ t~ s with 8 ss 0.15 for a fully occupied initial 
state. 

Following the arguments of in a related model, this 
observation can be used to provide an heuristic explana- 
tion for the fact that 8 ^ (3/v\\. 

It is known that, if the boundary acts as a sink or 
perfect trap (e.g. if p = 0), then, in a process starting 
with a fully occupied lattice, one observes a growing de- 
pletion zone around the boundary whose linear size l'(t) 
increases as I' ~ t° l , with a; = 1/2 (see [2(| and the 
next subsection). Thus, the density of active sites de- 
cays as t~ x / 2 . Hence, the influx of particles from the 
bulk to the leftmost site may be considered as an effec- 
tive time-dependent external field h(t) ~ t^ 1 / 2 . Making 
the assumption that this field varies so slowly that the 
response of the process (i.e. the actual average activ- 
ity at the boundary) behaves adiabatically, as if the field 
was constant, then in a critical process starting from an 
initially fully occupied state: 



po(t) ~ * ~ * 



-1/6 



(11) 



Owing to the time reversal property, this quantity should 
decay as the survival probability. This chain of heuristic 
arguments leads to the conjecture that the survival expo- 
nent is given by 8 = 1/6, in agreement with the numerical 
estimate 8 = 0.15(2). 

This unusual value of the exponent 8 is clearly related 
to the fact that the present problem is inhomogeneous. 
The argumentation presented above does not work for 
the CP, for example, since there is no special site and, 
therefore, a fully occupied lattice cannot be interpreted 
as a time dependent field acting on a special site. 
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G. Density profile 

Now, we consider the density profile p(x, t) in the bulk, 
where x £ N is the spatial coordinate (distance to the 
boundary), computed at the critical point. In the left 
panel of Fig.O we compare the data collapse of the curves 
p(x,t)t 1 / 2 as a function of x/t 1 / 2 with a Gaussian and 
observe an excellent agreement, indicating random- walk 
like behavior with a dynamical exponent z — 2. However, 
in contrast to a simple random walk, particles are mutu- 
ally correlated. This is illustrated in the right panel of 
Fig. O where the connected correlation function between 
two nearest neighbors: 

p pair (x,t) = (p(x+l,t)p(x,t))-(p(x+l,t))(p(x,t)) (12) 

in a system at the critical point is plotted against time. 
One observes an algebraic decay, a; -1 / 2 , with distance. 
According to the standard scaling theory this implies that 
f3/v± = 1/2 , confirming that z = v\\/vi_ = 2. Moreover, 
these results are in full agreement with field theoreti- 
cal calculations presented in Ref. [ijj (see section ITV C[) . 
which predict z = 2 and a — 1/2. 



III. MEAN FIELD APPROXIMATION 



Here, we study mean field approximations at differ- 
ent levels. Let us denote by rji the probability to find a 
particle at site i; the temporal evolution within a simple 
(one-site) mean field approximation is given by: 



dt 
dm 

dt 
drjj 

dt 



-0--p)vo + 7^1(1 - m), 



(13) 



p77o(l - T]l) + - (?72 + V0V1 ~ 2 Vl) i ( 14 ) 
2 (Vi+i + Vi-i - , for i = 2, 3, . . . (15) 



Note that the equations for the boundary site and its 
neighbor, Eq. IT3)) and Eq. (fT4|) . include quadratic terms 
due to the exclusion constraint, while the equation for 
sites at the bulk, Eq. (|15[) . describes in this approxima- 
tion a symmetric random walk, i.e. it is a diffusion equa- 
tion. The critical point within simple mean field theory 
(where the equation for 771 also becomes a diffusion equa- 
tion) is p c = 1/2. 

Considering a localized initial condition at the bound- 
ary, rji = 5ifi, after a transient time the densities at 
sites and 1 should, approximately, coincide. Therefore, 
from Eq. p3[) with r/o ~ f)i, it follows that, at criticality, 

In the stationary regime, Eq. (fT3"| leads to 7/0 ~ {p — 
1/2) for p > 1/2. From these results we have: 



, MF 



1/2. 



(3 



MF 



1 



(16) 



To obtain the survival exponent, 5, we follow the ar- 
guments of the preceding section and study the decay 



of activity from a fully occupied lattice, r\i = 1 for all 
i. Integrating Eqs. ((121), E q- 03) an d Ec l- C3 numeri- 
cally with this initial condition, we obtain an exponent 
in agreement with 



1/4. 



(17) 



A more accurate approximation can be obtained by keep- 
ing the correlation between the first two sites, which is 
expected to be more relevant than the correlation be- 
tween other neighboring sites. Such a pair- approximation 
was used recently in a model where a boundary site also 
plays a special role plj . In this approximation, the mas- 
ter equation reads: 



do 00 

dt 
dcrpi 

dt 
d(Jio 

dt 
do 11 

dt 

dr/2 

dt 

di]j 
dt 



(1 -p)<r w + ^[<7 i(i - m) - ^00772], (18) 

(1 -p)a n + 7^[oooV2 ~ CT oi(2 - m)}, 
-(7io + ^[o- i - owm + (Jii(l - 772)], 
po 10 - (1 - p)a n + i [cr 10 i]2 - cr n (1 - %)], 



(% + O-ll + ^01 - 2t7 2 ) , 

(rji+i + r)i-i - 2r)i) for i = 3,4, . . . 



where cr SQSl is the probability that the occupation num- 
bers of the first two sites are sq and s%. Numerical inte- 
gration of these equations leads to an improved critical 
point estimation, p c « 0.634, but to the same mean-field 
exponents as above. 



IV. RELATED MODELS AND FIELD 
THEORETICAL APPROACHES 



A. Bosonic variant 

The model defined above is fermionic in the sense that 
each site can be occupied by, at most, one particle. We 
now consider a bosonic variant without such a constraint. 
This means that diffusion is independent of the config- 
uration of particles and that particles can be created at 
the boundary site without restriction. More specifically, 
the update rules are: 

(a) A particle is chosen randomly. 

(b) If the particle is located at the leftmost site it can: 
create another particle at the leftmost site (sq = 
so + 1) at rate A, die (so = So — 1) at rate c, or 
diffuse to the next neighbor at rate D. 

(c) If the particle is located in the bulk, it diffuses to 
the right or to the left at equal rates D. 
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The corresponding master equation is: 
dP({n},t) 



dt 



= \[(n - l)P(n - 1, ...,*) - n P({n},t)] 
+ a[(n + l)P(n + 1, t) - n P({n}, *)] 
-l, nj + l,...,t) (19) 



D 



(y> 

+ P(...,ni + l, nj - 1, ...,t) - 2P({n},i) 



where P({n},i) is the probability to find a given con- 
figuration {n} = no,ni,n2 ■ ■ ■ and the sum runs over all 
nearest neighbors, j, of site i (recall that site has only 
one neighbor). Defining the state vector: 



|V(*)> = E p (t n >' t )K n »' 

{n} 



(20) 



where |{n}) = ®i\rii) denotes the usual configuration 
basis, the master equation can be expressed in the form 



dt 



\m) = -H\m), 



(21) 



where H is the time evolution operator. Using bosonic 
creation and annihilation operators, defined by di\rii) = 
n.i\n.i — 1) and S| = \fii + 1), the master equation 
Eq. (|2H)) can be shown to correspond to the time evolution 
operator: 



H = D^2{&\-&])(a i -& j ) 

m 

+cr(aj - l)a + AaJ(l - aJ)«o- 



(22) 



In this formalism, the expectation value of an operator 
B is given by (B) = (l|%(t)) where (1| = £{™L<MI- 
As is the case for the bosonic contact process [22|], the 
equations for the time evolution of the density of particles 
close. From the Heisenberg equation of motion, ^ = 
[H,B] and Eq. (|22]l. one obtains: 

^ = D( Pl po) + A Po (23) 
= D(p i+ i +pi_i - 2pi) i= 1,2,3... 

at 

where P i(t) — (al(t)ai(t)) = (a.i(t)} and A = A — a. 
Alternatively, one could have written a Langevin equa- 
tion equivalent to Eq. (|22p . and from it, averaging over 
the resulting noise, one readily arrives at the same set of 
equations Eq. (f2"3"l) . 

From these equations, we can see that the critical point 
is A = 0, where Eq. (|23"l) is a diffusion equation. In the 
continuum limit, Eq. (|23|) reads: 



dp{x,t) d 2 p(x,t) 



dt 



dx 2 



AS(x)p(x,t) 



(24) 



where x is the spatial coordinate and, without loss of 
generality, we have set D = 1. We note that in order 
to take the continuum limit in equation (|23p . a site — 1, 
with p-i = po, has to be introduced, so that appropriate 
boundary conditions are satisfied. The solution of this 
inhomogeneous diffusion equation is: 



p(x,t) 



5(QG(x,C,t)dC+ 



JO 



A5(() P (t,T)G(x,(,t~T),d(dT (25) 



where G(x,(,t) = ( e -(*+C) 2 /m + e -^~C) 2 /(4t))/(y^) 
is the Green function and the first term in the right 
hand side comes from the initial condition p{x, 0) = 6(x). 
From Eq. PS"]) we have 



wt 



2A 



dt- 1 / 2 



Po(t) 



(26) 



,— 1/2 

where po(t) — p(0, t), and the operator dt -i/ 2 ? defined by 



i/2 Jo y/w{t-T) 



dt 



(27) 



is a half integral operator [23[. Equation |26|) involves 
(owing to the delta function in the interaction term in 
Eq. [|24j0 only the density at the leftmost site. This jus- 
tifies the mapping of this model onto an effective one-site 
non-Markovian process (see next section). Using some 
rules for half integration [23| to solve Eq. ([26]) . we find: 



Po(t) 



4A exp(4A 2 t)erf(-2A%/i), 



(28) 



where erf(x) is the error function. This implies that, 
above the critical point, po grows exponentially in the 
long time limit, and does not reach a stationary value, 
i.e. there is a first order transition and, hence, /3 = in 
this bosonic model. From equation Eq. ([28]) . we deduce 



(3 1 = 1 and vn — 2. We have not been able to calculate 
the survival-probability exponent exactly, but numerical 
simulations suggest <5 = 1/4, in agreement with the mean 
field exponent. 



B. Partially bosonic variant 

Let us now introduce a partially bosonic variant of the 
previous model by retaining the exclusion constraint only 
at the boundary, but not in the bulk. The rules, in this 
case, are: 

(a) A particle is randomly chosen. 

(b) If it is at the leftmost site, it can generate a particle 
at site 1 (provided that s\ = 0) with probability p 
or die (sq := 0) with probability 1 — p. 
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(c) Particles in the bulk diffuse to the right or to the 
left with the same probability, 1/2. 

Numerical simulations show that this variant exhibits 
the same critical behavior as the original model, even 
if the critical point is shifted to p c = 0.6973(1). This 
shows that the fermionic constraint is relevant only at the 
boundary, where it induces a saturation of the particle 
density and leads the transition to become continuous. 



C. Models with pair annihilation at the boundary 

In the models discussed so far, particles at the bound- 
ary either create an offspring or die spontaneously at 
some rate. Instead, a very similar model was introduced 
in Ref. [l7j . for which particles at the boundary annihi- 
late only in pairs. In its fermionic variant, particles at 
sites and 1 annihilate with each other (provided that 
both sites are occupied) at some rate, while isolated par- 
ticles at the boundary cannot disappear: 

present models: A — > 2A , A — > , 
models of Ref. [17] : A -> 2 A , 2 A -> . 

Analogously, one can define a bosonic version, in which 
two particles at the boundary can annihilate. In the fol- 
lowing discussion we consider these two variants in d spa- 
tial dimensions where, as is the case d=l, only a single 
site has "special" dynamics. 

A detailed field theoretical analysis of these pair- 
annihilating models was presented in [l7| . In the bosonic 
case, proceeding as above (see Eq.(|22jl) one obtains the 
following time evolution operator: 

H = D^(at-a])(a i -a i ) 

(il) 

+<t[(4) 2 -l]ajj + AaJ(l-aJ)fio- (29) 

which, after eliminating higher order terms and taking 
the continuum limit, is equivalent to a Langevin equa- 
tion identical to the one for DP except for the fact that 
all terms, except for the Laplacian, are multiplied by a 
6 function at the boundary; i.e. the non-diffusive part 
of the dynamics operates only at the boundary. An e- 
expansion analysis of Eq. (j2"9"|) (see [13]) leads to a = 1/2 
and z = 2 as exact results in all orders of perturbation 
theory, and to (3 — 1 — 3(4 — 3d)/8, up to first order 
in e = 4/3 — d around the critical dimension d c = 4/3. 
Also, it was shown that the time reversal symmetry is 
preserved. 

We have verified all these predictions in computer sim- 
ulations of the bosonic annihilation model. For instance, 
from the time decay of po{t), as shown in Fig. [6l we 
determine S = 0.21(3), while from a finite size scaling 
analysis of the saturation values of the order parameter 
at criticality we measure 0/v± = 0.51(2) (see Fig. [5]), in 




FIG. 6: (Color online) Temporal behavior of po for the bosonic 
pair-annihilating model, starting from a homogeneous initial con- 
dition for different system sizes (from L = 64 to L = 2048). The 
exponent f3/v± can be measured from the scaling of the different 
saturation values as a function of system size (see inset; yellow 
line). Also, in the inset (dashed green line), we show the scaling 
of saturation values for the fermionic version of the same model, 
showing the same type of scaling. 



reasonable agreement with the expected results, 6=1/6 
and /3/vx = 1/2, respectively. Moreover, from spread- 
ing simulations (not shown) we estimate a ~ 1/2 and 
z 2. All the exponents are in agreement with the ones 
presented in the previous section for single particle anni- 
hilation models. 

Actually, a simple argument explains why the model of 
scctionlnland the pair-annihilation model share the same 
critical behavior. This is plausible because the chain re- 
action A — > 2A — > in the model with pair annihilation 
generates effectively the reaction A — > of the model 
considered with CP-like dynamics. 

Hence, the field theoretical predictions discussed above 
[HI, E3 a PPly a l so to the CP-like model. In d = 1, the 
one-loop prediction = 5/8 = 0.625 [l7j], is not far from 
the exponent measured in section [II] = 0.68(5). 

On the other hand, the fermionic version of the pair- 
annihilating model has been conjectured to yield in a 
different universality class, and a prediction for its criti- 
cal exponents is made in [17j (for instance, = 1). Our 
numerical simulations disprove such a claim; all the mea- 
sured critical exponents for the fermionic variant of the 
pair-annihilation model arc numerically indistinguishable 
from their bosonic counterparts (see Fig. [6]). 

In summary, all the defined models, either with single 
particle annihilation or with pair-annihilation, fermionic 
or bosonic, exhibit a boundary induced phase transitions 
and, except for one of them, they all are continuous and 
share the same critical behavior. The exception to this 
rule is the CP-like model without a fermionic constraint 
at the boundary, which lacks of a saturation mechanism 
in the active phase, leading to unbounded growth of par- 
ticle density at the leftmost site above the critical point 
and to a discontinuous transition. 
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V. RELATION TO A (0 + l)-DIMENSIONAL 
NON-MARKOVIAN PROCESS 

In Ref. [ItJ , by integrating out the fields related to dif- 
fusion in the bulk from the corresponding action, it was 
shown that the class of boundary-induced phase transi- 
tions into an absorbing state considered here can be re- 
lated to a non-Markovian single site process. The prop- 
erties of such a spreading process on a time line has been 
studied in further detail in Ref. [l8j |. 

On an heuristic basis, the relation can be explained as 
follows: consider the CP-like model only from the per- 
spective of the leftmost site. A particle at the origin may 
die or create a new particle that will go for a random walk 
coming back to the origin after a time r. What happens 
during this random walk is irrelevant from the perspec- 
tive of the leftmost site; the only relevant aspect is the 
time needed for a created particle to come back to the 
boundary. Once it returns it may die or create new off- 
springs which, on their turn, will undergo random walks 
in the bulk. 

Our simulations above show that the fermionic con- 
straint is irrelevant in the bulk. Therefore, we can con- 
sider without lost of generality the bulk-bosonic version 
in which there is no effective interaction among diffus- 
ing particles. In this case, the probability distribution 
of the returning time to the origin has the well-known 
asymptotic form (24[: 

P(t) ~ t- 3 ' 2 . (30) 

Taking all these elements into account we define the fol- 
lowing non-Markovian model on a single site [l6| : 

(a) Set initially s(t) := S tt o for all times, t. 

(b) Select the lowest t for which s(t) = 1. 

(c) With probability p, generate a waiting time r ac- 
cording to the distribution Eq. (|3"0")). truncate it to 
an integer, and set s(t + r) := 1; otherwise (with 
probability 1 — p) set s(t) := 0. 

(d) Go back to (b). 

The process runs until the system enters the absorbing 
state (s(t') = for all t' > t) or a predetermined max- 
imum time is exceeded. The density of particles at the 
leftmost site of the original model is related to (s(t)) in 
the single-site model, the survival probability at time t 
is given by the fraction of runs surviving at least up to 
t, and the initial condition s(t) := S t fi corresponds to 
start with a single particle at the boundary in the full 
model. Critical exponents can be defined as in the origi- 
nal model. However, the simulation results for the single- 
site non-Markovian model are more reliable because it is 
possible to perform much longer runs and, in the case 
of off-critical simulations, one can work with smaller val- 
ues of A. With time-dependent simulations at the crit- 
ical point p c = 0.574262(2), we obtained a = 0.500(5) 
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FIG. 7: (Color online) Off-critical data collapse with the one- 
site model: (s(t)) A -2 * 3 as a function of tA 413 for different 
values of A, with (5 = 0.71(2). 



and S — 0.165(3), in good agreement with the conjec- 
tured values a = 1/2 and 6=1/6. As an example, 
we show the results of supercritical simulations in Fig. 
where we obtained a convincing data collapse by plotting 
(s(t))A~ 2 @ as a function of tA 4 @ for different values of 
A with /3 = 0.71(2). The latter estimate is in agreement 
with j3 = 0.68(5), coming from the original model. 

As shown in previous studies (see e.g. [25[ and ref- 
erences therein), a non-Markovian time evolution with 
algebraically distributed waiting times P(t) ~ T -!- K jg 
generated by so-called fractional derivatives 9 t K which are 
defined by: 

A/||(«0 Jo 

where k € [0, 1] and A/ji (k) = — T(— k) is a normalization 
constant. Hence, we expect this model to be described by 
a DP-like 0-dimensional Langevin equation with a half- 
time derivative, instead of the usual one, to account for 
the non-Markovian character of the model: 

dl /2 p(t) = ap(t)-p(t) 2 +at) (32) 

where a is proportional to the distance from critical- 
ity and £ is a multiplicative noise with correlations 
(£(£)£(£')) = p(t)S(t — t'). This equation can be obtained 
from the effective action that arises when the fields re- 
lated to diffusion in the bulk are integrated out, and the 
relation of the order of the fractional derivative in a gen- 
eralized one-site model with the dimension in the full 
model is k — (2 — d)/2 [17]. An analysis of this one-site 
model with general k and a comparison with the results 
coming from field theory is presented in [l8j . 
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VI. CONCLUSION 

We have studied boundary-induced phase transitions 
into an absorbing state in one-dimensional systems with 
creation/annihilation dynamics at the boundary and sim- 
ple diffusive dynamics in the bulk. The non-trivial dy- 
namics at the boundary induces a phase transition in the 
bulk. We have analyzed such a transition for different 
though similar models, including different ingredients: ei- 
ther single-particle annihilation or pairwise annihilation, 
fermionic constraint or lack of it, etc. 

A particular bosonic version can be exactly solved; ow- 
ing to the lack of any saturation mechanism, the density 
of particles grows unboundedly in the active phase, lead- 
ing to a discontinuous transition with trivial critical ex- 
ponents. 

The rest of the analyzed models exhibit a continuous 
transition and define a unique universality class. At the 
bulk, the dynamics is governed by random-walks, entail- 
ing the exponent values z = 2 and a — 1/2. On the 
other hand, some critical exponents take non-trivial val- 
ues: i) the survival probability from a localized seed at 
the boundary exponent, which from an heuristic argu- 
ment supported by simulations results, turns out to be 
5 = 1/6, as well as ii) the order parameter exponent, 
f3 = 0.71(2). The remaining exponents can be obtained 
from these ones using scaling relations. 



Finally, it has been shown that the class of bound- 
ary induced phase transitions studied here can be related 
to a single-site non-Markovian process. This process is 
particularly suitable for numerical simulations and it is 
also of conceptual interest in the sense that it shows that 
nonequilibrium phase transitions can occur even in + 1 
dimensions by choosing an adequate non-Markovian dy- 
namics. It is also convenient for the comparison of the 
results obtained form the e-expansion and simulations 

mm- 

The models studied here possibly constitute the sim- 
plest universality class of nonequilibrium phase transition 
into an absorbing state, in the sense that the transition 
occurs because of the special dynamics of just one site 
and, in contrast to DP, some critical exponents can be 
obtained exactly from the field theory. 
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